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; Abstract 
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We calculate perturbative renormalization constants for the AB=2 four- 
^ ■ quark operators in lattice NRQCD. Continuum operators &7/x(l — 75)9 — 

, 75)5 and b{l — ^^)q 6(1 — 75)^, which are necessary in evaluating the mass 

^ \ and width differences in -B^^^^ — -B^^^^ systems, are matched at one-loop with 

CD ' corresponding lattice operators constructed from the NRQCD heavy quarks 

and the 0(a)-improved light quarks. Using these perturbative coefficients. 



o 

■ we also reanalyse our previous simulation results for the matrix elements 

, of the above operators. Our new results are free from the systematic error 

\ of 0{as/{aMb)) in contrast to the previous ones with matching coefficients 

evaluated in the static limit. 
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Typeset using REVTgX 



I. INTRODUCTION 



The B meson decay constant and the B parameter in the B — B mixing are crucial 
quantities for determining the Cabbibo-Kobayashi-Maskawa (CKM) mixing matrix elements 
\Vtd\ and \Vts\ from the experimental values of the oscillation frequency AMfi(s)- While the 
lattice calculation of the decay constant has reached a satisfactory level where the systematic 
error except for the quenching effect is about 10%, the B parameter Bb still has a large 
uncertainty of about 30% even in the quenched approximation Further effort in the 
lattice calculation is required to constrain the CKM matrix elements more tightly. 

In the limit of infinitely heavy quark mass, lattice calculation of the B parameter has been 
performed by several authors using the static action and the (9(a)-improved (or unimproved) 
light quark actions for which the perturbative matching factor of relevant four-quark 

operators in continuum and lattice definitions is available • "^^^ problem of large one- 
loop coefficient raised in Refs. |^,^ is not essential when used with the tadpole improved 
perturbation theory [jl2| as discussed in Refs. [^,|^ , where they find that the results of several 
groups are in reasonable agreement. 

The next step towards the final prediction is to incorporate the correction from finite 
b quark mass Mb, which can be systematically included using p/Mb expansion, where p is 
the typical momentum of the gluons and quarks. A naive order counting suggests that 
the correction is about AQco/Mb ~ 10%, when the size of the QCD scale is assumed to be 
around 350 MeV. On the lattice, the nonrelativistic QCD (NRQCD) ||13[ provides a necessary 
formulation to calculate the p/Mb corrections, and an exploratory lattice calculation of the 
B meson B parameter has already been made ||l^Jl5|] . One of the main drawbacks in that 



calculation is, however, that the one-loop coefficients for the infinitely heavy quark mass is 
used instead of those for lattice NRQCD. This approximation introduces a large systematic 
uncertainty of order as/{aMb), which is as large as 10-20% for a typical value of inverse 
lattice spacing 1/a ~ 2 GeV and almost equivalent to or even larger than the physical size 
of the p/Mb correction itself. 

Calculation of Bb using the relativistic lattice actions for heavy quark [16-23] is another 
possibility to study the finite heavy quark mass correction. These calculations, however, may 
suffer from large 0{aM) {0{asaM) or (9((aM)^) for the (9(a)-improved actions) systematic 
error and the uncertainty in the extrapolation to the b quark mass from lighter heavy quark 
masses, for which simulations are performed. 

In this paper, we compute the one-loop renormalization constants for AB=2 four-quark 
operators constructed with the NRQCD heavy quarks and with the (9(a)-improved light 
quarks on the lattice in order to remove the error of the order as/{aMb) in the lattice 
calculation of Bb- We consider the leading dimension six operators and neglect dimension 
seven operators which would remove errors of 0{asAQCD/Mb) or O^asaAqcD)- The one-loop 
coefficient for the dimension seven operators which corresponds to 0{asaAQCD) corrections 
has been obtained in Ref. in the infinitely heavy quark mass limit. 

Using the renormalization constants obtained in this work, we reanalyze the simulation 
data of Ref. |]14[ to obtain an improved result for Bb, which is free from the large systematic 
uncertainty of 0{as/{aMb)). The central value is increased by about 12% with this new 
analysis, which is within the size of errors expected by a naive order counting argument. 

Another important application of our perturbative work is the lattice calculation of 
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Bs, which is a S-parameter necessary to evaluate the width difference in the — B(^s) 
mixing. An exploratory lattice NRQCD study with the one-loop matching in the infinitely 
heavy quark mass limit is found in Ref. ||2^. We reanalyze the data in that work with the 
renormalization constant containing the finite heavy quark mass effect. We find that our 
new analysis resulted in a change of the value of the bag parameter, but it remains within 
the expected size of the error in the previous analysis as is also the case for Bb- 

This paper is organized as follows. We summarize the definition of the lattice NRQCD 
action in Sec. |I|, and the heavy-light four-quark operators in Sec. Sec. |^ is the main 
part of this paper, where we present the results of the one-loop matching calculations for 
bilinear operators ( [IV A[ ) and for the four-quark operators ( [IV B| ). The reanalysis of our 
previous NRQCD simulations are given in Sec. 0. Sec. ^ is devoted to our conclusion. 
Details of the one-loop calculations are collected in Appendices. The Feynman rules for the 
lattice NRQCD action is given in Appendix ^ while several expressions of one-loop integrals 
and amplitudes are summarized in Appendix ^ and y respectively. 



II. LATTICE FORMULATION OF NRQCD 

In this section we briefly summarize the definition of the NRQCD action used in the 
following perturbative calculations. A complete formulation of the lattice NRQCD is found 
in Ref. 

Our NRQCD action is defined by 



^NRQCD = - KQ){x,y)Q{y)+Y.xHx){l - K^){x,y)x{y). 



(2.1) 



x,y 



The nonrelativistic two-component spinor fields Q and x represent a heavy quark and an 
anti-quark respectively. Their evolution is described by [| 



Kqix.y) 

Kx{x,y) 



1 - 
1 



2n 



a5H" 



5i-)f/l 1 - 



5t^U, 1 - 




{x,y), 
{x,y), 



(2.2) 
(2.3) 



2n J \~ 2 y V" 2 / V" 2n 

where n denotes a stabilization parameter introduced in order to remove an instability arising 



from unphysical momentum modes in the evolution equation. Note that following Ref. |T3 



all the link variable f/^ in the NRQCD action is always divided by the mean field value Uq 
determined from the plaquette expectation value. This tadpole improvement will give rise 
to 0{g'^) counter term in the Feynman rule. 

The operator 5^^^^ is defined as 6^^\x,y) = 5x4±i,j/4'^x,y5 and the Hamiltonians Hq and 
6H are 



^The evolution equations ( [2.2[ ) and ( [2.3| ) are slightly different from the definition used, for example, 
in Ref. [p3], where the (1 — ai^o/2^^)" terms appear inside of the (1 — a6H/2) terms. 
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2aMo ' 

where aMg denotes a bare heavy quark mass in lattice unit. The operator A*^^) = Yl\=i ^f'^ 
is a Laplacian defined on the lattice through A • , the second symmetric covariant differen- 
tiation operator in the spatial direction i. The space-time indices x and y are implicit in 
these expressions. The Hamiltonian 5H represents the effect of the spin-(chromo)magnetic 
interaction, in which B is the chromomagnetic field defined as a standard clover-leaf oper- 
ator. (7 is a gauge coupling, and is a constant to parametrize the strength of the a ■ B 
interaction. It should be tuned until the NRQCD action reproduces the same dynamics as 
that of continuum relativistic action. We take the tree level value = 1. The relativistic 
four-component Dirac spinor field h is related to the two-component nonrelativistic field Q 
and X appearing in the NRQCD action in Eq. ( |2.1| ) via the Foldy-Wouthuysen-Tani (FWT) 
transformation 

Kx) = Y.R{^,^){^fll^. (2.6) 

where R is defined as 

R = l- ^ , (2.7) 
2aMo ' ^ ^ 

where Ap-* is the first symmetric covariant differentiaon operators in spatial direction. 

The Feynman rules derived from the NRQCD action in Eq. ( |2.1| ) and the FWT trans- 
formation (|2.6| ) are given in the Appendix 0. The light quark action is the (9(a)-improved 
Wilson action [^, and the gluon action is the standard plaquette action. The Feynman 
rules for light quarks and gluons are also summarized in the Appendix 0. 



III. OPERATORS 



The B parameters Bl^ and Bs are defined using the /S.B=2 four-quark operators h^^{l — 
75)'? ^7m(1 ~ 15)1 ci-nd 6(1 — 75)g 6(1 — 75)5 respectively. In the perturbative matching, 
however, we have to consider other operators which mix under the radiative correction. 
Since the lattice regularization violates the chiral symmetry, some operators that do not 
appear in the matching between continuum regularizations are also necessary. We define 
the following set of operators. 

OvLL = hixPiq h,^PLq, (3.1) 

OvRR = ht^PRQ hi^PRQ, (3-2) 



^We use a notation instead of the usual Bb in order to emphasize that it represents a matrix 
element of the "LL" operator. 



OvLR = h-i^Piq Iry^PRq, 
OsLL = bPLQ bPLQ, 
OsLR = bPLQ bPng, 
OvLL = hf^PiT^q b^^PLT'^q, 
OvRR = hf.PRT'^q h.PRT^q, 
OvLR = h^PiTy h^PRT'^q, 
OsLL = bPiT'^q IPlT^ 
OsLR = bPiT'^q bPnT'^q, 



(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 



where Pl and Pr are chirality projection operators Pl/r = (1 =F75)/2, and is a generator 
of the SU (N) group. The operators with a tilde contains a summation over the SU{N) gen- 
erators T". Fierz identities relate the 'tilde' operators in Eqs. (|3.6|) - (|3.1CI|) to those without 
tilde in Eqs. (O)-(O) as 



OvLL 
OvRR 
OvLR 
OsLL 
OsLR 



N 



2N 
N-1 



-O 



VLL, 



o 



VRRi 



o 



SLR, 



2N 

-l^OvLR 

T^J^^SLL - -^UvLL, 

-1^0 SLR - -^OvLR- 



(3.11) 
(3.12) 
(3.13) 
(3.14) 
(3.15) 



We use these relations to eliminate the 'tilde' operators from matching relations. We note 
that all equations except Eq. ( |3.14|) are exact, whereas Eq. (|3.14| ) is valid up to 0{p/Mq) 
correction terms described by dimension seven operators. When computing the matching of 
OvLL and Osll operators, the neglected terms give errors of C(asp/Mo) through one-loop 
mixing. 

In Sec.|V| we present our final result using the following set of operators in more conven- 
tional definitions 



Ol 


= &7m(1 - 75)g &7m(1 - 75)g, 




(3.16) 


Or 


= ^7m(i + i5)q ^7m(i + 75)g, 




(3.17) 


Os 


= -75)g &(i -75)g, 




(3.18) 


On 


= 267^(1 - 75)g 67^(1 + 75)g + 46(1 - 


75)g 6(1 + 75)g. 


(3.19) 


Om 


= 267^(1 - 75)g 67m(1 + 75)g - 46(1 - 


75)g6(l + 75)g. 


(3.20) 


Os 


= 6^(l-75)g^'&^(l-75)g\ 




(3.21) 


Op 


= 267^(1 - 75)g h^{l + -i^)q + ANhil 


-75)g 6(1 + 75)g. 


(3.22) 


Oq 


= 2^7^(1 - 75)g 67^(1 + 75)g + 46(1 


-l5)q b{^ + l5)q, 


(3.23) 


Ot 


= (2 + iV)67^(l - 75)g 67^(1 + 75)g - 


2(3iV2 -2N- 4)6(1 - 75)g &(1 + 75)?- 


(3.24) 



The indices i and j, which appear in the definition of O5, run over color of quarks, while 
other operators are products of color-singlet bilinear operators. As is obvious from Eqs. 
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( p.l| )-( pl24D , the set of operators in conventional definition are related to the first set of 
operators as 



Ol 


= 4 OvLL, 


(3.25) 


Or 


= 4 OvRR, 


(3.26) 


Os 


= 4 OsLL, 


(3.27) 


On 


= 8 {Ovlr + 20slr), 


(3.28) 


Om 


= 8 {OvLR — '^Oslr), 


(3.29) 


Os 


= 8 (^OsLL + 2iV^'^'^^) ' 


(3.30) 


Op 


= 8 {0vlr + 2N0slr), 


(3.31) 


Oq 


= 8iNOvLR + 2 0sLR), 


(3.32) 


Ot 


= 4(2 + N)OvLR - 8{3N^ -2N- A)Oslr. 


(3.33) 



IV. ONE-LOOP CALCULATION 



In order to match the operators defined in the continuum theory, say the MS scheme 
with the dimensional regularization, to the lattice counterparts, we compute the on-shell 
amplitude both in the continuum and on the lattice at one-loop level. 

Let Ox^{fi) and 0^*(l/a) be certain continuum and lattice operators defined at scale fi 
and 1/a respectively. The on-shell amplitude for a certain external state can be expressed 
by a linear combination of tree- level amplitudes (Oy)o, where the subscript Y runs over all 
possible operators which can mix with and O^* at one-loop, namely 



{Of{l/a)) = (Ox)o + ^EPxV(l/«)(Oy)o + 0(«.'), 



Y 



(4.1) 
(4.2) 



where = g'^/4:7i, and PxyM and px^yi^/^) represent the one-loop coefficients in the 
MS and the lattice schemes. We take zero spatial momentum on-shell free quarks for the 
external state. This choice is the easiest and sufficient to obtain the matching coefficients 
uniquely, since we restrict ourselves to the matching at lowest operator dimension, for which 
no derivative operator appears. 

Requiring that the both operators give identical one-loop on-shell amplitudes, we obtain 
the following matching relation 



Y 



0^^\l/a). 



(4.3) 



In the following we compute the coefficients Pxvif^) ^^"^ Px^yI^/^) heavy-light 
bilinear operators and AB = 2 four-quark operators. 
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A. Bilinear operators 



First of all, we give the expression for the matching of the bilinear operators for com- 
pleteness. Although the one-loop coefficients for the matching of the heavy-light vector 



and axial vector currents have already been obtained by Morningstar and Shigemitsu |28 



even through O^ap/Mo) and 0{aap), we present the one-loop matching coefficients for the 
general heavy-light bilinear operators for completeness. 

The one-loop expression of the perturbative on-shell amplitudes of the heavy-light bilin- 
ear operator hVq with arbitrary Dirac structure F is given as 



{{hVq){^)) 



(&rg)o, 



(4.4) 



for both continuum and lattice operators. There is no operator mixing in the lowest dimen- 
sion bilinear operators. In the continuum (the MS scheme with totally anti-commuting 75), 
the coefficient pp^'^(yu) is obtained |]7|, [T0| , [TT| as 



Ml 



2 3 3iJ2 

- - In — H HH' 



GH 



Mi 



11 

T 



(4.5) 



where Cp = [N'^ — 1)/2A^, and A denotes a gluon mass introduced to regularize the infrared 
divergence. The constants H, H' and G are defined through the following equations 



D JTT 

HV^Y.!,"^!,^ = — , GF = 74r74, 



(4.6) 



with space-time dimension D=4. The corresponding one-loop expression for the lattice 
operator is 



p'^\l/a) = G, 



|ln(aV) + l(G + C.) 



+{A'Kf[lA + GIb + {H- Gfic + {H- G){Id + If) + {H - G)GIe\ 



(4.7) 



The infrared divergence of form \Gf In is canceled between continuum and lattice expres- 
sions for any bilinear operator in the combination of the matching coefficient p^^^ — pp*. 
Numerical value of the light quark wave function renormalization factor Gi is 9.076 for the 
(9(a)-improved action. If one uses the normalization 1/ ^/u^ for the light quark field moti- 
vated by the tadpole improvement [|12|, the number becomes -0.164 for uq = (iTiUpy/^ 
(average plaquette), or 1.7106 for uq = l/Sncrit (critical hopping parameter). The heavy 
quark wave function renormalization Gh depends on the heavy quark mass oMq, and its 
numerical values are summarized in Table |. The constants I a, Ib, Ic, ^d, Ie and Ip are 
one-loop integrals from the vertex corrections shown in Figure ^ Their explicit expres- 
sions are given in the Appendix |B| (Eqs. (|Bl|) -( P6|) ), and their numerical values are given in 
Table |n|. 

In the following, we present the expressions of the matching factors for the temporal 
component of the axial current ^4 and for the pseudoscalar density P. 
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1. Axial vector current 



For the axial-vector current A4 with r=7574, we obtain H=2, H'=l, G=—l, for which 
the matching coefficients are 



PA4. - 



3, A2 3 

-- m — 

2 4 



(4.8) 



pf^il/a)=CF 

+ (47r)2(J^ - + 9Ic + 3Id - "He + 3/^ ) 



(4.9) 



does not have the logarithmic scale dependence because of the (partial-)conservation 
of the axial vector current. Combining the two expressions we obtain the matching relation 



A 



MS 



l + ^(21n(aX') + a)+0(«^ 



Af(l/a), 



where 



Ca — Cf 



I - \{Ch + Q) - {A7r)\lA -Ib + 9Ic + 3Id - 3Ie + 



(4.10) 



(4.11) 



Numerical values of the coefficient (a are listed in Table VII. 



2. Pseudoscalar density 

For the pseudo-scalar density P with r=75, we obtain if=— 4, H'- 
matching coefficients are 



3 In In — K + 



M2 



2~" Mi ' 4 



+ {4n)\lA -Ib + 9Ic - 3Id + 3/^^ - Sip) 



-1 and G=—l. The 



(4.12) 



Combining these expressions we obtain the matching relation 

,2 \ 



P""(/i) 



where 



13 1 

- - -(C, + Q) - {AnfilA -Ib + 9Ic - 3Id + 3/^; - 3/^) 



Cp — Cp 

Numerical values of the coefficient (p are listed in Table [VI| 



(4.13) 



(4.14) 



(4.15) 
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B. Four-quark operators 



Wc present the one-loop matching calculation of the four-quark operators Ovll and 
OsLL, which appear in the evaluation of the mass and width differences in the B^^s) ~ B^^s) 
systems. 



1. Ovll 

In the continuum theory preserving the chiral symmetry, the four-quark operator Ovll 
mixes with Osll under the radiative correction. At one- loop level, the on-shell amplitude 
of Ovll{.Ij) defined at scale is written as 



{O^Lii^)) = 1 + ^pTmVLLil^) 



MS 

47r' 



{Ovll)o + 



MS 

^^PVLL,SLL 



{Osll)o, 



where 



PyLW(/^) = 21n 



M-41n^ 
P' Mi 



35 



MS 
PVLL,SLL 



(4.16) 

(4.17) 
(4.18) 



The mixing coefficient Pvll,sll does not have the scale dependence at one-loop level. 

The same operator mixes with four operators Ovlr, Oslr, Osll and Ovrr on the lattice 
due to the lack of the chiral symmetry. We obtain the following expression for the on-shell 
amplitude of the lattice operator Ovll{^/(i)- 



where 



a 



'P^virVLLa/a) 



An 



{O 



VLL/O 



+ 

+ 



a 



« Jat 



An 



PvLLyLR 



a 



s Jat 



An 



PVLL,SLL 



{Ovlr)o + 

{Osll)o + 



a 



s Jat 



■^PVLL,SLR 



a 



s Jat 



An 



PVLL,VRR 



{Oslr)o 

{Ovrr)o, 



(4.19) 



PvLLyLLi^a) = -A\n{a'X') + Cp[Ci + C, 



PVLL.VLR 



Pv\l,slr = (47r)^ 



+ (47r) 
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I A + 2/c + AIe 



+ - (/g - 3(/^^ + // + /j + 21k) + 16/l + In) 



{An) 



2i^-2lB--{lD + lF) 
2{-AlB + YilD + lF) 



PVLL,SLL = W [-16(2/c - Ie)] , 



PVLL,VRR = (47r) 



4. 



(4.20) 

(4.21) 

(4.22) 
(4.23) 
(4.24) 
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The integrals Ia, Ib, Ic-, Id, Ie and Ip come from the diagrams in which a gluon mediates 
between heavy and hght quark hnes as shown in Figures |^ and These are the same 
integrals as in the vertex correction of the bilinear operators, whose numerical values are 
given in Table ||. Other integrals 7^, Ih, Ij, Ik, II, Im, In are characteristic of 
the corrections of the four-quark operators. The diagrams in which a gluon line mediates 
between two heavy quark lines (Figure ^ produce the five integrals Ig, Ih, h, Ij and Ik, 
whose expressions are given in the Appendix ^ (Eqs.(|B^)-( |Bll| )). Their heavy quark mass 
dependence is summarized in Table |T|. Other three, II, Im and 1^ defined in Eqs. (pi2|)- 



( [B14D correspond to the diagrams in which the two light quark lines are connected by a 
gluon line. These do not depend on the heavy quark mass, and their numerical values are 



-0.004635(3), / 



M 



-0.002433(1) 



'AT 



-0.012204(6). 



(4.25) 



In Appendix 0, one-loop expressions of the lattice on-shell amplitudes with general four- 
quark operators hVq hVq are presented. The above result in Eq. ( [4. 191 ) is obtained by 
applying the Fierz transformation for the color and spinor indices on the expressions 



Combining the continuum and the lattice results we obtain 



+ 
+ 



a 
Air 

Air 



O 



lat 
VLL 



lat 



PVLL,VLR 



l/a) + 



lat 



PVLL,SLR 



Ait 



{p^% 



SLL Pv\l,SLL 



An 

5LL(l/a) + 



l/a) 



Qlat 



a 



s Jat 



An 



PVLL,VRR 



/^lat I 
^VRR^ 



l/a). (4.26) 



The numerical values of Pvll,y listed in Tables |V. 

The matching relation for Ol is obtained using the conversion formula (|3.25|) - (|3.32|) as 
follows 



1 + 



a. 



An 



-21n^ + 41n(a2M2) + a,L 



+ ^CL,sOr{l/a) + ^CL,RO^^\l/a) 
+ ^CL,MO'^l{l/a). 



0't{l/a) 
CL,NO'^\l/a) 



(4.27) 



lat 



The coefficients Cl,s, Cl,r, Cl,n and (l,m are listed in Table [V^. The coefficient (l,m of 
vanishes in the static limit, and other coefficients agree with the previous work in the 

same limit. 



2. OsLL 

The matching relation for the operator Osll is obtained in a similar manner. 
The operator Osll mixes with Ovll with the radiative correction in the continuum. 
The on-shell amplitude with Osll for the zero momentum external state is written as 



{OITl{p)) 



'^s MS 



An 



PSLL,SLL ip) 



a 



{OsLl)o + -^PsLL,VLLip){OvLL)o, 



(4.28) 



10 



at one-loop level. The coefficients are 



MS , ^ 16, /i^ 4 



TFe 1/X 2 A 3 

The lattice operator O 



-In-^ + - In — ^ + -. 

3 Mi 3 2 

'^1^ mixes with five operators in the on-shell amplitude 

{Osll)o 



(4.29) 
(4.30) 



as 



1 + ^PsUsLLi^a 



+ 
+ 



-^PSLL,VLR {OvLr)o + 



a 



'^P''SLL,SLR 



{Oslr)o 



where 



-^P^SLiyRR {Ovrr)o, 

pfLL,sLLiM^) = 'l Ha'X') + C4Q + C] 

+ (47r)2 ^{AIa + 52Ic + 20Ie 

-2{Ig + Ih + Ii + Ij + 2Ik + In)) , 
PsLyLLiya) = hnia^X') + (An)' ^(2(-3J^ - 7Ic + Ie) 
-3Ig + Ih + Ii + Ij + 2Ik - IQh - 3/jv 



PsLsLR = W 

p'sLyLR = (4^)' 
PsLyRR = (4vr)^ 



^llB + filn + lF)] 
2[\ib - ^{Id + If)] 



-I 



M 



Numerical values of p''sll y given in Table 0. 

Combining the continuum and the lattice results we obtain 



+ 
+ 
+ 



1 + ^ {p¥LL,SLLip) - PsLL,SLLi^/' 

PsLyLL^/^ 



{ps^LLyiL^p) 



OfL^^la) 



lat 

^^PSLL,SLR 

(^s lat 
'^PsLLyLR 



OTlr^/^) 



(4.31) 



(4.32) 



(4.33) 
(4.34) 
(4.35) 
(4.36) 



follows 



+ [-^PslLyRR] O^iLyRRi^a). (4.37) 

The matching relation for O5 is obtained using the conversion formula ( |3.25|) - (|3.32|) as 

ows 
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+-- 

4:71 



^ l^ln-^ + - Inia^M^) + Cs s 
1, /i^ 2 



Or(l/a) 



(4.38) 



The coefficients ^5,5, Cs,l, Cs,r, Cs,p and Cs,t are listed in Table [VII1| . The coefficient Cs,t of 
O^* vanishes in the static limit, and other coefficients agree with the previous work p^j2l 
in the same limit. 



V. PHYSICS RESULTS 

Using the one-loop coefficients obtained in this work, we reanalyze the lattice NRQCD 
calculations of Bl [Q and of Bs ||2^. These previous calculations were performed with the 
lattice NRQCD action for heavy quark, but the perturbative matching of the four-quark 
operators were done using the coefficients in the infinitely heavy quark mass limit. Due to 
this approximation for the matching coefficients, the previous results contain errors of order 
as/{aM), which is one of the largest uncertainties among all the systematic errors. 



A. Bl 

The B parameter -Bl is defined through 

BM . , (5.1) 

where the scale fi is usually set at the b quark mass Mh- In the following analysis we use /i 
= 4.8 GeV. On the lattice we measure the parameters' 

U/aj - |(50|^iai(i/^)|o)(o|Af (l/a)|fiO)' ^ ^ 

for four-quark operators Ox=Ol, Os, Or, On and Om using the NRQCD action in Eq. 
( p. 11) . We performed the simulations on a quenched 16^x48 lattice at j3=5.9. Other details 
of the lattice calculations are found in Ref. [|1^ . 



The perturbative matching relation for the continuum operator Ol in Eq. ( [4.27|) may 
be used to obtain 

BM = ZL,LMfi,a)B'^\l/a) + Zl,s/a^B's"\'^/(^) 

+ZL,R/A^B'^\l/a) + ZL,N/A^B'^\l/a) + Zl^m/a^B'^}{1/ a), (5.3) 

where the matching factors are 
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ZL,L/A^^, a) = 1 + ^ (^2 In ^ + Cl,l - 2Ca) , (5.4) 

Zl,s/a^ = ^Cl,5, (5.5) 

(y. 

Zl,r/a^ = ^Cl./j, (5.6) 

The one-loop coefficients (l,x are defined in Eq. ([4.27| ) and plotted in Figure |^ as a function 
of l/aMg. The coefficient of the leading contribution (l^l — '^Ca becomes larger in magnitude 
toward lighter heavy quark. The mass dependence of Cl,l—'^Ca is relatively smaller than that 
of Cl,l itself, due to the cancellation of singlet diagrams (Figure |^) against the contribution 
of the denominator —2(a. Two mixing coefficients Cl,5 and Cl.tv become smaller when 1/M 
correction is incorporated. It is also important that Cl,m, which vanishes in the static limit, 
becomes non-zero for finite heavy quark mass. 

The matrix elements on the lattice i?^*(l/a) measured in Ref. |0 are shown in Figure |^ 



as a function of inverse meson mass 1/Mp. Their mass dependence is qualitatively well 
described by the vacuum saturation approximation. 

For the coupling constant as we choose the l^-scheme coupling ay(g*) with q* = 
2/ a. To estimate the size of higher order perturbative errors we also analyze with q* = 1/a 
and IT /a. 

Combining Z^ x/A'^ and S^* we obtain the contribution of each term Z l^x/a'^ B''-^^ {1 / a) in 
Eq. (^.3|) as shown in Figure p. In spite of the large mass dependence of the coefficients (l,s 
and Cl,n, there is no significant mass dependence in the corresponding combined quantities 
ZL,x/A^B'--^^{l/a), since those are canceled by the large and opposite mass dependence in 
B^g^ll/a) and i?j^*(l/a). Small increase of Zi^^m/a^B^^}{1/ a) from zero in the static limit is 
observed, which reflects the increasing trend of both Z^^m/a^ and B^^}{l/a). 

Total result for BL^rrih) is presented in Figure ^ by filled circles. Because of the cancella- 
tion of the large mass dependences in x and in i?^*(l/a), there is little 1/Mp dependence 
in our final result (ffiled circles). A small increase toward larger 1/Mp comes from the con- 
tribution of Oj\/. In this plot our estimate of systematic uncertainty is shown by error bars. 
The horizontal ticks attached to the error bars represent the size of statistical error, which 
is much smaller than the systematic errors especially for large 1/Mp points. 

We also plot our previous analysis with the same NRQCD action but using the pertur- 
bative matching in the static limit (open circles in Figure There is a small negative 
slope in 1/Mp so that the previous result is about 12% smaller than our new result at the B 
meson mass. An estimation of 0{as/{aMb)) error in our previous analysis is around 10%, 
when we assume an order counting argument with typical value of the strong coupling con- 
stant as ~ 0.3. In addition to this error, there are also other errors of 0{al), O^o^Aq^jj) 
and O^asaAqcD)- Thus the shift of our result does not exceed the systematic uncertainty 
discussed in Ref. fl^ . 

Our final numerical result is 

^^(mfc) = 0.85 ±0.03 ±0.11, (5.9) 
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where the first error is statistical and the second is systematic. The systematic error is esti- 
mated using the order counting of neglected contributions. In our new analysis in which the 
(9(a<j/(aMb)) error is removed, the remaining sources of uncertainty are the dicretization 
errors C(a^AQ^j^) (~ 5%) and O^asaKqcD) (~ 5%), as well as higher order perturba- 
tive error 0{a1) (~ 10%). Higher order contributions in the nonrelativistic expansion are 
OipisJ^QCD/Mb) (~ 2%) and 0{Aq^^/M^) (< 1%), which are not dominant uncertainties. 
We assume as ~ 0.3 and Aqcd ~ 350 MeV when we estimate the errors listed above, which 
are added in quadrature to give the systematic error of about 13% in the final result. 

The result (|5.9|) may be compared with the recent lattice calculations with relativistic 



heavy quark actions: 0.92(4)1;] and 0.93(8)lij;[] where the first error is statistical 
and the second is their estimate of systematic errors. It is encouraging that our result 
agrees with these relativistic calculations within the large systematic uncertainty in { f).9\ j. 
Although the systematic error in the relativistic results seem much smaller, it should be 
noted that the quoted systematic uncertainty could be underestimated. They extrapolate 
their simulation results performed around charm mass regime assuming the 1/M scaling 
without considering 0((aMo)^) errors. However, the 1/M dependence of the simulation 
results could be distorted by the C((aMo)^) error, which can be as large as 30% toward the 
heavier side in the naive order counting. In order to have a reliable prediction of Bi{mb), one 
has to at least include 0((aMo)^ error when extrapolating in 1/M, or take careful contiuum 
limit before doing 1/M extrapolation. Furthermore, the heavy quark expansion becomes 
questionable to describe the heavy quark in the charm quark mass regime when truncated 
at 1/M or at 1/M^. Therefore, an alternative method to fit the results would be to include 
the result in the static limit in order to constrain at least the leading term in the /IM 
expansion. 

Finally, we also obtain chiral breaking effect on the ratio of as 

^1^4— = 1-01 ± 0-01 ± 0-03, (5.10) 
for which the relativistic results are 0.98(3) [^] and 0.98(5) |23[| . 

B. Bs 

The B parameter Bs is defined as 

f(50|PA/S(^)|0)(0|pM5(^)|50) ^ ' 

We also define a ratio of the matrix elements of bilinear operators 

{0\Af^\B^) 



7^(^) 

and calculate 



(0|PA^^(/x)|50) 
Bsif^) (P0|0p(^)|i?0) 



(5.12) 



(5.13) 
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which is necessary in evaluating the Bg meson width difference . 
In the lattice simulation we measure 



for four quark operators Ox = Os, Ol, Or, Op and Ot- Note that the denominator of B 



I lat 
X 



is different from Eq. (|5.11|). The perturbative matching relation for the continuum operator 



O^^ in Eq. ([4.38|) may be used to obtain 



+Zs,R,A^B'^''\l/a) + Zs,P/A^B'p"'\l/a) + Zs^T/A^B'T^^'il/ a), (5.15) 



where 



^s,sM2(/i,a) 


= 1 + ? 


3 


Zs,L/Ai (/^, o) 


= ^filn 


Mi 










as . 




Zs,T/A'^ 


as . 





Mi 

-h 

3 



,2 i\/r2\ 



,2 T,j2\ 



:s,s-2uy (5.16) 

(5.17) 

(5.18) 
(5.19) 
(5.20) 

The one-loop coefficients defined in Eq. ( |4.38| ) are plotted in Figure |10| as a function of 
1/aMo. The I/Mq dependence is quite large for the leading contribution Cs,s — '^Ca, because 
the denominator in Eq. ( [5.13[ ) is not a vacuum saturation of the numerator and the cancel- 



lation of color singlet diagrams (Figure |^) does not take place. Among other coefficients, 
C,s,L has relatively large I/Mq dependence and the mixing of operator Ol becomes smaller 
as 1/aMo increases. 

The matrix elements i?^*(l/a) are shown in Figure [ll|. The 1/Mp dependence in Bs, 
Bp and B-p is significant, which is well described by the vacuum saturation approximation 
as discussed in Ref. ||2^. The contribution of each term Zsx/A'^B''-^^{l/a) in Eq. ( |5.15|) to 



B s{fTib) / Tl{mi,Y is plotted in Figure |T^, in which no clear 1/Mp dependence is observed. 

Total result for Bsimh) /TZ{mi,)'^ is presented in Figure |13| by filled circles. As is evident 
from the plot of each term Zs,x/A^B''-^^{l/a) (Figure |12|), there is no clear trend in the 1/Mp 
dependence of Bs{mb) /TZ{mbf' . Our previous analysis [0] with matching coefficients in the 



static limit is plotted with open symbols. Reduction of the result with the correction is quite 
large (~ 20%), but consistent with our estimate for the collection of systematic errors of 
C(a,/(aM)), 0{al), O^o^K^cd)^ O^aaKQco) (~ 20%). Main effect comes from the large 
1 / aMo dependence of C,s,s — 2Ca shown in Figure |10[ 
Our final numerical result is 

= 1.24 ± 0.03 ± 0.16, (5.21) 
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where the first error is statistical one, and the second is our estimate of systematic uncer- 
tainty obtained as in the analysis of B^. For the width difference we obtain 



0.107 ± 0.026 ± 0.014 ± 0.017 



(5.22) 



using Eq.(9) of Ref. [^. Errors are from the Bg meson decay constant Jb^ = 245(30) MeV, 
which is taken from the current world average of unquenched lattice calculations [|I|, from 
Bs{fnb) /TZ{'mb)'^ , and from an estimate of higher order contribution in the l/rrib expansion 



VI. CONCLUSIONS 

In this paper we have performed one-loop calculations of matching coefficients for AB=2 
four-quark operators defined using lattice NRQCD. This calculation allows to remove one 
of the dominant systematic errors characterized by 0{as/{aMb)) from the lattice simulation 
of the B parameters B^ and Bs- We find sizable I/qMq dependence in several one-loop 
coefficients, which affects the mass dependence of the B parameters as well as their absolute 
values at the b quark mass. 

We have also presented a reanalysis of our previous simulations and obtained results for 
Bl and for Bs/TZ'^ with reduced systematic error. The difference from our previous results is 
consistent with the estimate obtained with order counting argument. Remaining systematic 
uncertainty is dominated by unknown two-loop matching coefficients. 
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APPENDIX A: LATTICE NRQCD FEYNMAN RULES 



In order to simply the expression, we set the lattice spacing a = 1 throughout Appendix A 
and B. When deriving the Feynman rules from the NRQCD action, we followed the method 
which is explained in Ref. We also note that the Feynman rules for 0(1/M) NRQCD 
action with slightly different definition from ours are given in Ref. 



1. Functions 



We define the following functions which appear in the Feynman rules below. 

P^E4sin^| (Al) 
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= 1 



1 ^ /• 



nMo — i 



f^,M) = sin COS y 
JuAQi^ Q2) = COS 1 — - — I sm ( — — | sm 



(A2) 

(A3) 
(A4) 

(A5) 



1 [cos (l)^cos (g, + + cos (g^ + |)^cos (|^^ 
+ cos {q, + I) cos (g, + I) - cos (|) cos (| 



(A6) 



We also define 

v4f^ (/',/; Pi, Ail, P2,/i2) 

A^^^(/',/,Pl,/il,P2,/^2) 



2raMn 



«-i 



-1 



Li=o 

n-l 



4nMo 
1 



^^(0)(;'y^(0)(^)n-l- 



,i=0 























cos(^ 2 





Ml 



Ml 



(A7) 
(A8) 



i=0 j=0 



n—2—i- 



(2nMo)2 



M2 



'4Mn 



dHP{l'J;pi,iJ.i,p2,H2) = -i 
where S* denotes a four-by-four matrix 

= 



E^ ^'J>i /jti (^'i ) ' 
3 

jj=i 
3 



4M 



i=i 



(A9) 
(AlO) 

(All) 

(A12) 



a' 
a* 



Using these functions the Fourier components of the evolution operator in Eq. 
written as 

,(1) 



(A13) 
13) is 



+^(0) (;')"^(o) ^;y"c<(/', , ^^)] 5^^^ 

-^e-^('^+'^)AW(0"A(°)(/)"54,„ 



(A14) 
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+ e-*('-^^)M(i) (/', / - P2; pi, - p2, 1; P2, fi2) 

+ {l'rC{l', I - P2)A^'^ (/ - P2, 1; P2, fi2)dHi'^ (/', I - P2; Pi, 1^1) 
+ A^''\IYC{1 - p2, l)A^'\l', I - P2;pi, i^i)dHi'\l - p2, l;p2, 1^2) 
+ A^'^ {irC{l', l)dHP (/', l;pi,i^i,p2, 1^2) 

+ e-'(^~P^^^A^°\irA^'\irdHi'\l\l-p2;pi,i^i)dHi'\l-p2j;p2,^^^ 
+ [(e-''^ (/')" + e-'M(o) {1)^) A?^ (/', pi, //i, p2, //2) 



+^(0) (;/)n^(0) (;)n^^^.^ ^)^^(2) /■ /^2) 



^4/ii5/ii/i2) 



(A15) 



where 54^^ = 1 — 5, 



■4m- 



2. Feynman Rules 

We summarize the Feynman rules used in our calculation. 

In our convention the arrows in the heavy (anti-)quark propagator represent the flow of 
momentum irrespective of whether it is particle Q or anti-particle %. Other notations are 

• a, b, color index of quarks, 

• a, P, spin index of quarks, 

• /, /': momentum of quarks, 

• A, B, color index of gluons, 

• n, 1/, spin index of gluons, 

• k, ki, momentum of gluons. 

Double lines denote the heavy (anti-) quark propagators. 

We also need Feynman rules for heavy-light bilinear and four-quark operators. In this 
appendix, we give the Feynman rules for AB = — 1 heavy- light current as an illustration. 
Feynman rules for other operators can easily be deduced. 

In the diagrams involving a heavy-light current, heavy quark is incoming into the current, 
and heavy anti-quark is outgoing from the current. Light quark represented by a single line 
is, on the other hand, always outgoing. 

• Gluon propagator: 

^K^SSJ)SI^J)J)J)J)S^ ^^^.^(A;) = (5^%,^^^ (A16) 
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• Light quark propagator: 



a, a 



Jab 



a/3 



where S{1)-^ = ^ sin^ + (U' 



(A17) 



• 0{g) vertex for hght quark: 



a,a 




-t-f^ cos 



— r sm 



^ A=l J a/3 



^ (A18) 



• 0{g^) vertex for hght quark: 

The vertex from the clover term does not give any contribution to the diagrams we 
compute, thus we do not give the explcit expression here. 



a, a 




-|-{T^T^}a. 



r cos 



'I' + V 



Ml 



a/3 



+contribution from the clover term 



(A19) 



• Heavy quark propagator: 



a, a 



Jab 



1+74 

2 /a/3 
1+74 



Q{l) 



>ah 



2 /a/3 l-e-*'4y4(o)(0 



2n 



(A20) 



0{g'^) counter term introduced for the tadpole improvement: 

This term appears because we devide all the link variables 11^ in the NRQCD action 
by the mean field value of mq = 1 — g^u^'^^ . 



a,a 



Mo '"'I 2 y.^ 



(A21) 
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0{g) vertex for heavy quark: 



a, a 




v^^\l',l;k,i,) 



1+74 



a/3 



0{g'^) vertex for vertex: 




a,a 



' ((T^T^),,4')(/',/;A;i,/xi,A;2,/i2) 
+ {T''T\,v'^^\l', I; k2, //2, h, /xi)) - 



Heavy anti-quark propagator: 



Sab[^-^] Q{1) 



a, a 



2 J a/3 
1-74A Sab 



2 ya/3 l-e"*'4^(o)(0 



2n 



C(gr) vertex for heavy anti-quark: 




heavy-hght current: 



,1+74" 



a/3 
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• heavy-light current with the FWT rotation: 



-* HLx^ 



a, a 



2Mr ^ ' 



i=l 



Jab 



(A27) 



a/3 



vertex from rotated heavy-hght current: 



(1, A 
^ k 



a,a fo,p 



(, k\ 1 + 74 

1 , ^ 7u cos / 

271^0^^^ V 2 



{T\,S,, (A28) 



a/3 



(anti-) heavy- Ught current: 



a, a 



r I 



,1-74 



Jab 



a/3 



(A29) 



(anti-)heavy-hght current with the FWT rotation: 



a, a 







r ^ . , 1 - 74 



Jab 



r I 



a/3 



(A30) 



vertex from rotated (anti-) heavy-hght current: 



H, A 
^ k 

I 



a, a 



1 7u cos t H — 

2Mo''' V 2)^ 2 



{TXbK. (A31) 



a/3 



APPENDIX B: ONE-LOOP INTEGRALS 



We hst the integrals appearing in the one-loop calculations with the NRQCD action. 



G{l)S{~l)Q{~l) 



X 



(27r)4 

cos ^ sin k + sin I ) + ( sin^ + ? sm' '-^ ] ^X(/) 



i=l 



h] 1 

2 / 2' 
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I 



l + -^n E^in'^^E cos 



2 h 



<i=l 



i^sin' h cos' ^-^]Y{1) 

i=l 



dH 9(1 -P) 



(27r)4 (^2)2 ' 

"^'^ G{l)S{-l)Q{-l) 



(27r) 



X 



1 -r) I A C. 



—i — sin — sin I4 + -/ cos — 



+i f sin^ — ^ sin^ /j ] sin /4X(/) 

^ G(05(-0Q(-0 



3 



3 



^ 1 + ^/^ E-^/^E 



G(/)5(-/)g(-/) 



i=i j=i 
1 



2 ^7 

cos — 
2 



Esin^iiCOS^IlnO 

i=l 



6M0 



X 



E sin^ ^- sin ^ + ^ cos ^ sin Z{1) 



+ E sin' k [sin' | - ^ sin' /4 ) X{1) 



- ^ - Csi„ E "^os' 7^ E sin' - E ^i"^ ^' 



2 

(2^ 



- G{l)S{-l)Q{-l) 



i=l 



6Mn 



1 + %^^[' 
4 



X 



cos^Z'(Z) — ^sinZ4X(/)^ sin' Ij 



cos 



2 '' t 



+ i I E sin' h E cos' 7^ — E sin' cos' 7^ 1 sin U Y{1) 
,i=i j=i 2 2 

rf''/ 1 
G{l)S{-l) 



(27r) 

xEf^~ 



12Mn 



1=1 



2 2 ^« • 2 7 I • 2 7 2 • 2 7 
cos — — sm li + Csw sm / j cos — — Cgw 2^ sm 

2 2 



cos 



2 



2 ^ 



/ 



^G(/)gH)' (Esin^^X(/f + 



+ 



f dH 


2Mo 


'^(1-H 


1 (27r)4 


L-2iMoZ4 + ^'J 


^2 



sinh ^ 



(47r) 



2Mo' 
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Ih = 
h 



^ \i=l .7=1 



i=l j=l 



— sin^ Ij 



3^3 ^.3 
/j E cos^ 77 ~ E ^^^^ ~ E ^^^^ 

2 .^-^ 2.^^ 



X 



•3 i=i \i=i i=i ^ i=i / . 



1 / 1 1 

X — > sm cos — — > sm In, cos — 
12 I 2 ^ 2 

^ J (27r)4 4(f-j^ 

4 ; 1 4 



r dH e{l-P) 
^ J (27r)4 4(/2)2 ' 

1 r 4 



1 

^4 



- E sin^ If, 



2 1 2 • 2 7 • : 

- T^sw E E sin _ 

;=1 u=l p=l 



1 ~ ^ 



/„ sin' „ 
".=1 2 4 

i^/,E«in^^.cos^^^-i(P)^E' 

4 ' 1 ~ ^ I 

-c^j'^ ^sin^ I ^ cos^ , 



2 7 2 



| + >E«in^^, 
/i=i 



' cos^ ■ 

2 



! 7 2 

i/=l ^ 



^(/YE^ 



1~ ^ 1 

--/^Esin^L-^ 
, 2 ^4 

4 4/4 7 

2 7 A • 2 7 2 ''M 

2^ cos — - 4 2^ sm COS — 



■ 2 ''M 

sm — 



1 / * 
- sin 



where functions -^(/), and Z(Z) in the integrands are 
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n-1 



X{1) 



Z{1) = 



e^'M(°)(/)" + l 

Hi _i_ 1 



iC4 

4iWn' 
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(B15) 

(B16) 

(BIT) 
(B18) 



There are infra-red divergences in the integrals I a ( |B1| ), Iq ( p7|) and //, ( pi2|) , for which we 
subtract an continuum expression from their integrand in the region < 1. We, then, add 
back their analytic integral except for the In(aA) term, so that Ix becomes finite. When 
those integrals appear in the expressions of on-shell amplitude, the infra-red divergences will 
be added. 



APPENDIX C: ONE-LOOP FOUR-QUARK AMPLITUDES 



The lattice one-loop expression of the perturbative on-shell amplitudes is 



{bTq hVq) 



^lat ^lat 



{bTq bTq)o 



Q^s / -.^singlet 
^yj- \' ^ heavy— light 



I 'v'octet I -yoctet _i_ -yoctet 

heavy— light ~^ heavy— heavy ~^ light— light) i 



(CI) 



where Z\°-^ and Z]^^ are light and heavy quark wave function renormalizations respectively, 
and give by 



7lat 



zlat 



ln(a"A") + Ci 
-2ln(a^X^) +Ch 



(C2) 
(C3) 



The vertex corrections X's are classified by the topology of Feynman diagrams. Figure 
^ shows the diagrams in which the gluon line connects heavy and light quarks and the flow 
of color is closed. The amplitude of these diagrams is denoted as ^ heavy-light- Figure |^ 
the gluon line connects heavy and light quarks, but the color flow is not closed, which we 
call -^heavy-iigM- Figurcs ^ and | represent the diagrams in which the gluon line mediates 
between two heavy quarks or between two light quarks respectively. The color flow cannot 
close in these diagrams, and we denote them as Xf^f^^y-heavy 
The expressions of the one-loop amplitudes are the following. 



and Xf^M-jig/^t, respectively. 



^ heavy— light 



2 1 



1 



167r2 
3 _ 

+2^c J2 {Hlj'^ljliQ bTq)o 



\n{a'X') {bTq bTq)^ 
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i=l 



•y octet 

heavy— light 



+2/E^(6747ir7j74g hTq)o 
1 



(47r)2 

+2/B(6r74r«g 674rr"g)o 



2 ( /a - ^ ln(a2A2) ) {hVTy hVT'^q)^ 



+2Ic E (^r7,7,T'^g lrfaj'^T^q)o 
+2(/c + /f) E(^r7,T'^g 67irT"g)o 



octet 

heavy— heavy 



i=l 

+2/s^(6r7,74T'^g 674%rT»g)o 
j=i 

(47r)^ (Jg - 2^ ln(a2A^)) (6rT»g 6rT»g)c 

3 _ 

1=1 



"V^ octet 
light— light 



+ Ij + 2Ik) Y.{hV^'T-q hV^'T'^q)^ 
1=1 

(47r)2 (/L + ^ln(a2A2)) ^^(6r7^7,r«g 6r7^7, 

+/ME(6r7^T'^Or7^T'^g)o 
+lN(bTT''q 6rT"g)c 
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TABLES 



aMo 


n 


Ch 


12.0 


2 


2.94(12) 


10.0 


2 


2.63( 7) 


7.0 


2 


1.86( 5) 


6.5 


2 


1.67(12) 


5.0 


2 


0.89(13) 


4.5 


2 


0.52(12) 


4.0 


2 


0.07(11) 


3.8 


2 


-0.13(12) 


3.5 


2 


-0.49(11) 


3.0 


2 


-1.20( 9) 


2.6 


2 


-1.93( 8) 


2.1 


3 


-2.97( 8) 


1.5 


3 


-5.10( 6) 


1.3 


3 


-6.10( 6) 


1.2 


3 


-6.67( 6) 


0.9 


4 


-8.68( 5) 



TABLE I. Wave function renormalization constants for the NRQCD heavy quark 



aMo 


n 


/a 


Ib 


Ic 


Id 


Ie 


If 


12.0 


2 


0.030932(15) 


-0.016562(9) 


0.000794(1) 


-0.000029(0) 


0.001020(1) 


-0.000519(0) 


10.0 


2 


0.030277(15) 


-0.016041(8) 


0.000936(1) 


-0.000031(0) 


0.001209(1) 


-0.000623(0) 


7.0 


2 


0.028754(14) 


-0.014812(8) 


0.001282(1) 


-0.000032(0) 


0.001678(2) 


-0.000890(0) 


6.5 


2 


0.028388(14) 


-0.014512(8) 


0.001367(1) 


-0.000031(0) 


0.001791(2) 


-0.000958(0) 


5.0 


2 


0.026993(13) 


-0.013325(7) 


0.001702(2) 


-0.000024(0) 


0.002256(2) 


-0.001246(0) 


4.5 


2 


0.026371(13) 


-0.012802(7) 


0.001853(2) 


-0.000019(0) 


0.002469(2) 


-0.001384(0) 


4.0 


2 


0.025697(24) 


-0.012179(7) 


0.002033(2) 


-0.000010(0) 


0.002723(3) 


-0.001557(0) 


3.8 


2 


0.025327(12) 


-0.011894(6) 


0.002116(2) 


-0.000006(0) 


0.002843(3) 


-0.001639(0) 


3.5 


2 


0.024791(12) 


-0.011422(6) 


0.002252(2) 


0.000003(0) 


0.003039(3) 


-0.001780(0) 


3.0 


2 


0.023692(12) 


-0.010492(5) 


0.002524(2) 


0.000025(0) 


0.003423(3) 


-0.002076(0) 


2.6 


2 


0.022733(23) 


-0.009579(5) 


0.002793(3) 


0.000052(0) 


0.003826(4) 


-0.002396(0) 


2.1 


3 


0.021122(11) 


-0.008532(5) 


0.003253(3) 


0.000111(0) 


0.004524(4) 


-0.002966(1) 


1.5 


3 


0.018483(9) 


-0.006337(3) 


0.003995(4) 


0.000252(0) 


0.005707(6) 


-0.004153(1) 


1.3 


3 


0.017439(9) 


-0.005419(3) 


0.004320(4) 


0.000331(0) 


0.006230(5) 


-0.004792(1) 


1.2 


3 


0.016842(17) 


-0.004921(3) 


0.004502(4) 


0.000381(0) 


0.006529(6) 


-0.005191(1) 


0.9 


4 


0.014482(17) 


-0.003598(2) 


0.005228(5) 


0.000609(0) 


0.007825(7) 


-0.006921(1) 



TABLE II. Integrals Ia-, Ib-, Ic, Id-, Ie and Ip- 
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^ i\/r 
alViQ 


n 


r 


T 




T 

h 




12. U 


2 


A Aiocro/o\ 

0.01858(8) 


A AAAA'7r' / r\\ 

0.000075(0) 


A AAAAOA/a\ 

0.000029(0) 


A AAA1/^A/a\ 

-0.000160(0) 


A AAAAAO/'a\ 

0.000008(0) 


lU.U 


2 


0.01674(7) 


0.000104(0) 


0.000041(0) 


-0.000231(0) 


0.000013(0) 


7.L) 


2 


0.01z23(5j 


0.000193(1) 


0.000073(0) 


-0.000471(0) 


0.000037(0) 


D.5 


2 


A A1 "tA'7/r'\ 

0.01107(5) 


A AAAOAA/ 1 \ 

0.000220(1) 


A AAAA01 f r\\ 

0.000081(0) 


A AAAr'/1'7/A\ 

-0.000547(0) 


A AAAA/ir/'A\ 

0.000045(0) 


5.U 


2 


0.00doo(3) 


A AAAOOO/l\ 

0.000338(1) 


A AAA1 1 O / A\ 

0.000118(0) 


A AAAAO A I A^ 

-0.000924(0) 


A aaaaac:/'a\ 

0.000095(0) 


4.5 


2 


A AA/10A/o\ 

0.00439(2) 


A AAA/1AA/o\ 

0.000400(2) 


A AAA1 O /I / 1 \ 

0.000134(1) 


A AA1 1 AC\(C\\ 

-0.001140(0) 


A AAA10'7/'a\ 

0.000127(0) 


4.0 


2 


A AA1 00/0\ 

0.00183(2) 


0.000474(0) 


A AAA1 Cro/A\ 

0.000153(0) 


A AA1 A A^ { C\\ 

-0.001443(0) 


A AAA1'7/^/a\ 

0.000176(0) 


o o 
O.O 


2 


0.0U0d3(z) 


0.000515(2) 


0.000160(1) 


-0.001599(0) 


0.000203(0) 


d.5 


2 


A AA1 OO / 1 \ 

-0.00138(1) 


A AAArOA/'o\ 

0.000580(2) 


A AAA1 '71 /" 1 \ 

0.000171(1) 


A AA1 OOCr/A\ 

-0.001885(0) 


A AAAOC A /A\ 

0.000254(0) 


d.U 


2 


-0.00545(2) 


A AAA*?! 0/'0\ 

0.000718(3) 


A AAA"! Or/'l \ 

0.000185(1) 


-0.002566(1) 


A AAAOOr/^A\ 

0.000385(0) 


2.6 


2 


-0 00964(2) 


000862(1) 


000182(0) 


-0 003416(1) 


000563(0) 


2.1 


3 


-0.01528(5) 


0.001154(5) 


0.000146(1) 


-0.005237(1) 


0.000977(1) 


1.5 


3 


-0.02721(9) 


0.001680(5) 


-0.000234(1) 


-0.010264(3) 


0.002270(1) 


1.3 


3 


-0.03274(9) 


0.001932(8) 


-0.000613(2) 


-0.013665(4) 


0.003201(2) 


1.2 


3 


-0.03592(4) 


0.002081(2) 


-0.000926(2) 


-0.016038(4) 


0.003865(2) 


0.9 


4 


-0.04662(12) 


0.002760(9) 


-0.002752(9) 


-0.028512(8) 


0.007507(4) 


TABLE III. Integrals Iq, Ih, Ii, Ij and Ik for Cs^=l. 


aMo 


n 


J,at 

PVLLyLL 


PVLLyLR 


Jiat 

PVLL,SLR 


PVLL,SLL 


^lat 

PVLL,VRR 


12.0 


2 


19.78(16) 


10.751(2) 


20.353(4) 


-1.4375(4) 


0.51217(9) 


10.0 


2 


19.10(9) 


10.480(2) 


19.588(4) 


-1.6789(4) 


0.51217(9) 


7.0 


2 


17.47(7) 


9.840(2) 


17.740(4) 


-2.2483(5) 


0.51217(9) 


6.5 


2 


17.07(16) 


9.687(2) 


17.292(4) 


-2.3830(5) 


0.51217(9) 


5.0 


2 


15.47(17) 


9.091(2) 


15.511(3) 


-2.9032(6) 


0.51217(9) 


4.5 


2 


14.74(16) 


8.828(2) 


14.705(3) 


-3.1306(7) 


0.51217(9) 


4.0 


2 


13.85(15) 


8.519(2) 


13.742(3) 


-3.3948(7) 


0.51217(9) 


3.8 


2 


13.46(16) 


8.381(2) 


13.301(3) 


-3.5138(8) 


0.51217(9) 


3.5 


2 


12.78(15) 


8.153(2) 


12.568(3) 


-3.7104(8) 


0.51217(9) 


3.0 


2 


11.44(12) 


7.710(2) 


11.103(3) 


-4.0868(9) 


0.51217(9) 


2.6 


2 


10.10(11) 


7.286(2) 


9.641(2) 


-4.4438(9) 


0.51217(9) 


2.1 


3 


8.26(11) 


6.894(1) 


7.778(2) 


-4.9772(10) 


0.51217(9) 


1.5 


3 


4.76(8) 


6.059(2) 


3.902(1) 


-5.7882(12) 


0.51217(9) 


1.3 


3 


3.27(8) 


5.773(2) 


2.154(2) 


-6.0962(13) 


0.51217(9) 


1.2 


3 


2.47(8) 


5.642(2) 


1.155(1) 


-6.2536(13) 


0.51217(9) 


0.9 


4 


0.06(7) 


5.596(1) 


-2.098(1) 


-6.6608(15) 


0.51217(9) 



TABLE IV. Matching coefficients Oyi^y for Csw=^- The infrared divergence — 41n(a^A^) is 
subtracted from Py^j^ VLL- The tadpole improvement is applied such that the normalization of 
light quark field becomes ^/^Kcrit- 
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^ 7\/f 

alVI 


n 


PSLL.SLL 


PSLLyLL 


PSLL,SLR 


^lat 

PSLL.VLR 


PSLL,VRR 


"ion 

12. U 


2 


15.26(10] 


-1.8396(16] 


-8.3367(14] 


-1.2359(3] 


-0.12804(3] 


inn 
lU.U 


o 
z 


15.5U(yj 


1 Toon/ 1 o\ 

-1.7380(13] 


-8.1851(13] 


1 1011 /'o\ 

-l.loii(2j 


-0.12804(3] 


T n 
7.U 


z 


16.U9(7j 


-1.4903(9] 


-7.8368(15] 




n 1 oon /I /OA 
-0.12804(3] 


0.5 


2 


16.24(16] 


-1.4299(8] 


-7.7565(15] 




n 1 oon /I /OA 
-0.12804(3] 


K n 


2 


"1 O "I / "I 

16.81(17] 


-1.1879(6] 


^7 /I cnn / 1 /I \ 

-7.4509(14] 


-0.8857(2] 


-0.12804(3] 


4.5 


2 


17.1)5(16] 


-1.0782(5] 


T oon /I / 1 \ 

-7.3204(13] 


-0.8267(2] 


-0.12804(3] 


A n 
4.U 


2 


1 / .35(15 ] 


-0.9466(5 ] 


- / .1 ( 14(lo ] 


-0. (55/(2] 


-0.12804(0 ] 


o o 
O.O 


2 


17.5U(16] 


-0.8862(4] 


T 1 n^ji /I OA 
-7.1061(13] 


-0.7230(2] 


n 1 oon /I /OA 
-0.12804(3] 


o c 
O.O 


2 


17.71(15] 


n TO ceo/ A \ 

-0.7852(4] 


^7 nm / "I o\ 

-7.0013(12] 


n alio A /OA 
-0.6684(2] 


n 1 oon/1 /OA 
-0.12804(3] 


o.U 


z 


18.15(12] 


-0.5837(3] 


-6.8068(11] 


-0.5587(2] 


-0.12804(3] 


2.6 


2 


18.59(11) 


-0.3823(3) 


-6.6347(11) 


-0.4481(1) 


-0.12804(3) 


2.1 


3 


19.56(11) 


-0.1036(3) 


-6.5966(10) 


-0.2982(1) 


-0.12804(3) 


1.5 


3 


20.93(8) 


0.4334(2) 


-6.4940(10) 


0.0127(1) 


-0.12804(3) 


1.3 


3 


21.57(8) 


0.6692(3) 


-6.5593(10) 


0.1589(1) 


-0.12804(3) 


1.2 


3 


21.96(8) 


0.8006(3) 


-6.6362(10) 


0.2443(1) 


-0.12804(3) 


0.9 


4 


23.95(7) 


1.2199(4) 


-7.3536(9) 


0.5465(1) 


-0.12804(3) 



TABLE V. Matching coefficients O^'^l^ y for Csw=^- The infrared divergences — | In(a^A^) and 
I In(a^A^) are subtracted from p^^^, gj^j^ and (^slj^ vll respectively. The tadpole improvement is 
applied such that the normalization of light quark field becomes y/8Kcrit- 



aMo 


n 




Cln 


Clm 


Cls 


Clr 


12.0 


2 


-31.46(16) 


-5.232(3) 


-0.1441(0) 


-6.572(2) 


-0.51217(9) 


10.0 


2 


-30.78(9) 


-5.068(3) 


-0.1720(1) 


-6.326(2) 


-0.51217(9) 


7.0 


2 


-29.14(7) 


-4.677(2) 


-0.2426(1) 


-5.755(2) 


-0.51217(9) 


6.5 


2 


-28.75(16) 


-4.583(2) 


-0.2604(1) 


-5.623(2) 


-0.51217(9) 


5.0 


2 


-27.14(17) 


-4.211(2) 


-0.3343(1) 


-5.105(2) 


-0.51217(9) 


4.5 


2 


-26.41(16) 


-4.045(2) 


-0.3693(1) 


-4.878(1) 


-0.51217(9) 


4.0 


2 


-25.53(15) 


-3.847(2) 


-0.4126(1) 


-4.613(2) 


-0.51217(9) 


3.8 


2 


-25.13(16) 


-3.757(2) 


-0.4330(1) 


-4.494(2) 


-0.51217(9) 


3.5 


2 


-24.46(15) 


-3.609(2) 


-0.4676(1) 


-4.299(2) 


-0.51217(9) 


3.0 


2 


-23.12(12) 


-3.315(2) 


-0.5400(2) 


-3.923(2) 


-0.51217(9) 


2.6 


2 


-21.78(11) 


-3.026(2) 


-0.6169(2) 


-3.564(2) 


-0.51217(9) 


2.1 


3 


-19.93(11) 


-2.695(2) 


-0.7515(2) 


-3.029(2) 


-0.51217(9) 


1.5 


3 


-16.45(8) 


-2.002(1) 


-1.0267(3) 


-2.221(2) 


-0.51217(9) 


1.3 


3 


-14.95(8) 


-1.712(1) 


-1.1739(3) 


-1.913(2) 


-0.51217(9) 


1.2 


3 


-14.15(8) 


-1.555(1) 


-1.2660(3) 


-1.756(2) 


-0.51217(9) 


0.9 


4 


-11.74(7) 


-1.136(1) 


-1.6612(4) 


-1.351(3) 


-0.51217(9) 



TABLE VI. Matching coefficients C for Ovll 
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aJViQ 






n 




Ca 




Cp 


12. U 






2 




-14.59(8) 




1 1 O/^ / o\ 

-11.26(8) 


lU.U 






2 




-14.23(5) 




1 1 0^ / C \ 

-11.26(5) 


7.L) 






2 




1 o ooio\ 

-13.33(3j 




-11.29(3) 


6.5 






2 




-13.11(8) 




-11.30(8) 


5.U 






2 




1001/ n\ 

-12.21(9) 




11 OA/ r\\ 

-11.34(9) 


4.5 






2 




1 1 'VA / 0\ 

-11.79(8) 




-11.36(8) 


A A 

4.0 






2 




-11.28(7) 




1 1 Ar\/''7\ 

-11.40(7) 


o o 
O.O 






o 
z 




-11.05(o) 




-11.40(o) 


d.5 






z 




-10.66(7) 




-11.42(7) 


O A 






z 




-9.86(6) 




11 A r / /^\ 

-11.45(6) 


2.6 






2 




-9.04(5) 




-11.51(5) 


2.1 






3 




-7.89(5) 




-11.89(5) 


1.5 






3 




-5.47(4) 




-12.27(4) 


1.3 






3 




-4.31(4) 




-12.49(4) 


1.2 






3 




-3.63(4) 




-12.67(4) 


0.9 






4 




-1.14(3) 




-13.66(3) 






TABLE VII. Matching 


coefficients Ca,p 


for A4,P 






n 








Csp 


CST 


CSR 


1 O A 

Iz.U 


2 


-5.26(16) 




3.3396(16) 


A ^Cy<AC/lA\ 

0.65405(10) 


A A1 A A t A / 1 \ 

-0.014414(1) 


A 1 00 A /I /0\ 

0.12804(3) 


1 A A 

lU.U 


2 


-5.50(9) 




3.2380(13) 


0.63362(10) 


A A1 'TO A ^ / I \ 

-0.017204(1) 


A 1 00 A /0\ 

0.12804(3) 


^ A 

7. (J 


2 


-6.09(7) 




2.9903(9) 


0.58467(9) 


-0.024260(2) 


0.12804(3) 


6.5 


2 


-6.24(16) 




2.9299(8) 


0.57291(9) 


-0.026043(2) 


A 100A/1/'0\ 

0.12804(3) 


C A 

5.U 


2 


-6.81(17) 




2.6879(6) 


0.52643(8) 


-0.033428(2) 


A 100A/1/'0\ 

0.12804(3) 


4.5 


o 
z 


-7.05(16) 




2.5782(5) 


0.50564(oj 


-0.036930(2j 


0.12804(o) 


4.U 


o 
z 


-7.35(15) 




2.4466(5) 


0.48095(7) 


-0.041261(3) 


A 1 OOA /I /'0^ 

0.12804(3) 


o o 
O.O 


2 


-7.50(16) 




2.3862(4) 


0.46972(7 j 


-0.043297(3j 


0.12804(3) 


6.0 


2 


-7.71(15) 




2.2852(4) 


A /I CT 1 1 /I /" 

0.45114(7) 


-0.046764(3) 


A 100A/1/0\ 

0.12804(3) 


O A 


2 


-8.16(12) 




2.0837(3) 


0.41443(6) 


A AT/IAAI /o\ 

-0.054001(3) 


0.12804(3) 


2.6 


2 


-8.59(11) 




1.8823(3) 


0.37830(6) 


-0.061688(4) 


0.12804(3) 


2.1 


3 


-9.56(11) 




1.6036(3) 


0.33695(5) 


-0.075151(5) 


0.12804(3) 


1.5 


3 


-10.93(8) 




1.0666(2) 


0.25032(4) 


-0.102676(7) 


0.12804(3) 


1.3 


3 


-11.57(8) 




0.8308(3) 


0.21406(4) 


-0.117394(8) 


0.12804(3) 


1.2 


3 


-11.96(8) 




0.6994(3) 


0.19437(3) 


-0.126601(9) 


0.12804(3) 


0.9 


4 


-13.95(7) 




0.2801(4) 


0.14209(3) 


-0.166128(11) 


0.12804(3) 



TABLE VIII. Matching coefficients C for Osll 
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1/aM^ 

FIG. 6. l/aMo dependence of the one-loop coefficient for the matching of Ol- For Ql we plot 
C,l,l — '^Ca (cicle), and others are C,l^s (diamond), C,l,r (triangle right), Cl.a^ (triangle up) and Cl,m 
(triangle left). 
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FIG. 7. Matrix elements B^^\l/a) in Eq. ( |5.2D measm'ed on the lattice. Dashed curves repre- 
sent a quadratic fit to the data, and extrapolated static limit is shown by a open symbol. 
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1/Mp (GeV ^) 

FIG. 8. Contribution of individual operators Z ^ x/A'^ / to Bl{ii). Symbols are X = L 
(circle), 5 (diamond), R (triangle right), N (triangle up), M (triangle left). 




1/Mp (GeV ^) 

FIG. 9. 1/Mp dependence of Bi{rrn,) obtained with the NRQCD action (filled circles). The 
same calculation but with the matching coefficient in the infinite mass limit is shown by open circles. 
Other symbols in the static limit are calculations with the static action by UKQCD |^ (triangle 
up), Gimenez-Martinelli Q (triangle right) and Christensen et al. Q (square). We reanalysed their 
raw data with the same method as ours, namely we used Eqs. ( |5.3D -( ^^ with ai/(2/a) and C's in 
the static limit. 
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1/aM^ 

FIG. 10. 1/aMo dependence of the one-loop coefficient for the matching of Os- For C,s we plot 
C,s,s — 2Ca (cicle), and others are C,s,L (diamond), Cs,r (triangle right), C,s,P (triangle left) and C,s,T 
(triangle down). 



2.0 r 
1.6- 
^ 1.2f- 



0.8 



0.4 
0.0 

Or 

-4 ■ 



-12 [ 

53 

Q, 

«3 -16 



0.0 



-20 
-24 



0.0 



o 



o 



o 



0.1 

1/Mp (GeV~') 



0.2 



O 



O 



o 



o 



o 



0.1 

1/Mp (GeV~') 



0.2 




l/Mp(GeV ) 



FIG. 11. Matrix elements i?^*(l/a) in Eq. ( ^.141) measured on the lattice. 
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FIG. 12. Contribution of individual operators Zgx/A'^^x^i^/^^) to Bgip) /TZ{^j)'^ ■ Symbols are 
X = S (circle), L (diamond), R (triangle right), P (triangle left), T (triangle down). 
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FIG. 13. 1/Mp dependence of Bs{mi,) /n{mbY obtained with the NRQCD action (filled cir- 
cles). The same calculation but with the matching coefficient in the infinite mass limit is shown 
by open circles. 
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